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We show how to construct fully symmetric, gapped states without topological order on a honey¬ 
comb lattice for S = 1/2 spins using the language of projected entangled pair states (PEPS). An 
explicit example is given for the virtual bond dimension D = 4. Four distinct classes differing by 
lattice quantum numbers are found by applying the systematic classification scheme introduced by 
two of the authors [S. Jiang and Y. Ran, Phys. Rev. B 92, 104414 (2015)]. Lack of topological 
degeneracy or other conventional forms of symmetry breaking, and the existence of energy gap in 
the proposed wave functions, are checked by numerical calculations of the entanglement entropy 
and various correlation functions. Our work provides the first explicit realization of a featureless 
quantum insulator for spin-1/2 particles on a honeycomb lattice. 


Introduction - A modern theme of much interest in 
condensed matter systems is the classification of possi¬ 
ble phases of quantum matter in low dimensions. First 
noted in the context of quantum Hall physics, it has be¬ 
come clear that different quantum phases are labeled 
often by their topological characters rather than bro¬ 
ken symmetries as in the conventional Ginzburg-Landau 
paradigm [H. How to define such quantum orders and 
classify states accordingly in a precise way has intrigued 
theorists for several decades. 

A powerful guide in the classification effort is the “no- 
go” theorem such as the celebrated Lieb-Schultz-Mattis 
theorem in one dimension Q and its higher-dimensional 
generalizations due, for instance, to Oshikawa Q and 
Hastings [3| , stating that lattice spin models having a 
half-integer spin per unit cell must remain gapless, or if 
gapped, would either break conventional symmetries or 
turn into a topological state with fractionalized excita¬ 
tions. Powerful as they are, though, integer spin systems 
are not covered by these theorems. In one dimension we 
have some well-established results for integer-spin chains, 
e.g. S = 1 Haldane spin chain, saying that the ground 
state can be both gapped and featureless. 

Turning to two dimensions, search for an analogous 
featureless phase not addressed by the no-go theorems 
can best proceed by an explicit identification of *5 = 1 
models on a Bravais lattice or S = 1/2 models on a 
honeycomb lattice where an even number of sites form a 
unit cell [5f0. One such construction was given recently 
with S = 1 model on a square lattice [H , while attempts 
to construct featureless states for S' = 1/2 spins on a 
honeycomb lattice has met with partial success so far ^ 
0 ]. In this paper, we provide an explicit construction of 
the spin-1/2 wave function on a honeycomb lattice that 
preserves the full set of lattice symmetries plus time- 
reversal and SU(2) spin rotation, in addition to being 
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FIG. 1. (Color online) Schematic figures of the graphical rep¬ 
resentation of (a) tensor network on honeycomb lattice (b) 
u- and v-site tensors comprised of three gray legs and up¬ 
ward black leg describing Hilbert spaces of three virtual spins 
and a single physical one respectively and bond tensor (green 
ellipses) connecting neighboring virtual spins. Bravais unit 
vectors are chosen as ai = x and a 2 = (x + V^y)/2. 


devoid of topological order. We dub such a state feature¬ 
less quantum insulator, or FQI. Lacking both symmetry 
breaking and topological order, such states do not per¬ 
mit a straightforward field-theoretic description [HI . We 
instead use the recently developed classification scheme 
of the tensor network wave functions |8l 4l5|. in particular 
the one proposed by two of the authors [8|, to identify 
all possible spin-1/2 FQIs on the honeycomb lattice for 
a given bond dimension of the tensor network. Inten¬ 
sive numerical check carried out by the authors confirm 
that the proposed FQI state is indeed devoid of any con¬ 
ventional order, and has the topological entanglement 
entropy [iEIil of zero. 

Spin-^ symmetric PEPS on honeycomb lattice - The 
honeycomb lattice we work on and various notations for 
site and bond labels are shown schematically in Fig.lTJa). 
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To construct a honeycomb PEPS, we associate every 
site/bond of the honeycomb lattice with a site/bond ten¬ 
sor. A site tensor is formed by one physical leg which 
supports physical spin-1/2 degrees of freedom, and three 
virtual legs, while a bond tensor is formed by two vir¬ 
tual legs forming a nearest-neighbor bond as shown in 
Fig.IHb). Since every leg is associated with a specific lo¬ 
cal Hilbert space of spins, a tensor can be viewed as a 
quantum state in the Hilbert space of the tensor product 
of all its leg Hilbert spaces. The physical wave function 
is obtained by contracting all connected virtual legs of 
site tensors and bond tensors. 

In order to construct a fully symmetric, topologically 
trivial, and energetically gapped state with spin-1/2 per 
site on a honeycomb lattice, we rely on the recently devel¬ 
oped classification algorithm of quantum phases in terms 
of PEPS proposed by two of the present authors in Ref. 
0. Let us briefly review the procedure here. Due to the 
existence of an enlarged Hilbert space by virtual legs, the 
mapping from site/bond tensors to physical wave func¬ 
tions are actually many to one. Namely, a wave function 
which is globally symmetry-preserving can have its con¬ 
stituent site/bond tensors “gauge-transformed”, by act¬ 
ing with an arbitrary invertible matrix V{s^i) on a virtual 
leg i of the site tensor at s and simultaneously with its 
inverse matrix on the corresponding leg of the 

bond tensor. There may also exist special gauge trans¬ 
formations that leave every site/bond tensor invariant, 
up to U(l) phase factor. Those transformations form a 
group, named the invariant gauge group (IGG), which 
governs the low ener gy g auge dynamics of the state as 
shown in Refs.RSl andliol. 

Let us consider various symmetry operations on PEPS. 
According to the general remarks above, invariance of the 
physical wave function j?/) (up to a U(l) phase) under a 
specific symmetry operation g implies the following gen¬ 
eral transformation rules for the site and bond tensors: 


= e,W, {go T^) 

T^ = Wg{goT^). ( 1 ) 

We label the site tensor at s as , and as the tensor at 
the bond b. Wg is a leg-dependent gauge transformation 
acting on virtual legs of tensors. Here go can be spin ro¬ 
tation, time reversal or any lattice symmetry operations 
on a physical Hilbert space. Symmetry implementation is 
projective in the sense that the operation on the physical 
indices by g can be “compensated for” by the gauge op¬ 
erations Wg on virtual indices and the U(l) phase factor 
0^. The site-dependent phase factor 0^ also allows us 
to c^ture the symmetry quantum numbers of the state 

Eollowing the framework developed in Ref.[8|, symme¬ 
try group operations pertinent to the particular lattice 
geometry can be cast as a set of algebraic equations. By 


solving them, one obtains highly constrained forms of 
all the gauge transformation matrices Wg and 0^ associ¬ 
ated with the physical symmetric operation g. We should 
mention that IGG will in general enter these algebraic 
equations, influencing the outcome of the solutions for 
Wg and 0^. In keeping with the spirit of the present pa¬ 
per, which is the search for EQI in the case of the spin-1/2 
honeycomb lattice, we set IGG to be trivial, i.e. as an 
identity element. 

Details of the classification procedure for symmetric 
PEPS with trivial IGG and how to solve for the kE^’s 
and are found in the Supplementary Information 

(SI) [2^. We should emphasize that the final expression 
for the EQI tensors are quite transparent and can be un¬ 
derstood without the full knowledge of the classification 
scheme. In the end, we obtain 

WTAs,i) = WTAs,i) = \ 

WcAu,a/b) = WcA'^A) = I, W^Cs(m,c) = 

(u,a) = W^, (u, h/c) = xc, W,, 

Wa{v, a) = XCe XaCe Wa, Wa{v, b/c) = Wa, 

M 

Wr{s, i) = WT = 0(ld, ® 

k=l 

M 

WenisA) = Wen = 0(ld. ® 

k=l 

®Ti(<s) = 0T2(«5) = I, QCeiu) = XCg^ 

= 1, = XCeXaCe: ©T = 1- (2) 

Some words on the notations are in order. Two sublat¬ 
tice sites within the unit cell (yellow dashed box in Eig.[T]) 
are labeled s = u, v, while i = c label the three vir¬ 
tual legs associated with each lattice site as indicated in 
Eig.[T](b). Various symmetry operators considered here 
are translations by lattice vectors ai and a 2 , denoted 
respectively as Ti and T 2 , six-fold rotations about the 
center of the hexagon (Ce), mirror reflection a along the 
dashed line in Eig.[Tl time-reversal symmetry (T), and 
spin rotation symmetry. Each virtual leg has the Hilbert 
space consisting of M different species of spins, each la¬ 
beled as 5/c, I < /c < M. Eor each spin Sk one further in¬ 
troduces the “flavor” degeneracy of for a total virtual 
Hilbert space dimension D = <^/c(25'/c + !)• To have 

trivial IGG, we are required to assign only half-integer 
spins at the virtual legs [ 2 ^ . We further have 

M _ 

VF<, = 0(w^®l2S.+i), (3) 

k=l 

where W^ is a d/^-dimensional real matrix satisfying 
{W^Y — ^dk' Notice that all kk^’s in Eq. ([2]) are trans- 
lationally invariant (independent of s). 
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( XCs, XctCs ) 

(+ 1 , + 1 ) 

(+ 1 ? ~ 1 ) 

(“ 1 ^ + 1 ) 

(- 1 ,- 1 ) 

Ce 

+ 1 

+ 1 

-1 

-1 

(7 

+ 1 

-1 

-1 

+ 1 


TABLE L Lattice quantum numbers {xcq , XctCq) for Cq and a 
operations for four different classes of featureless states. Sign 
of the wave function |'0) changes by the amount shown in the 
second and the third rows under the Ce and a operations, re¬ 
spectively, for each state characterized by the pair of quantum 
numbers (xce,XcrCe) in the first row. 

According to our analysis 0 , there are four differ¬ 
ent symmetric featureless PEPS classes characterized by 
XCe = and XaCe = ^^5 regardless of the bond dimen¬ 
sion D. The PEPS wave functions belonging to different 
classes indeed can be distinguished by their lattice quan¬ 
tum numbers obtained for the torus geometry with odd 
number of unit cells, as shown in TablelH For instance, a 
PEPS wave function with (xCe^XaCe) — (+ 1 ? in¬ 

variant under the Cq operation and gains an extra phase 
— 1 under the a operation. After these formal matters, 
the remaining task is the explicit construction of site and 
bond tensors and the examination of physical properties 
for the state obtained from contracting the site/bond ten¬ 
sors. It should be cautioned that, even when the PEPS 
wave function is seemingly invariant under all symmetry 
operations, there is a chance that it actually describes 
a spontaneous symmetry breaking phase. To rule out 
these possibilities and to ensure that the constructed 
PEPS state is indeed a symmetric FQI, one should care¬ 
fully measure the correlation functions for varying system 
sizes. 

We will focus on a particular case where every vir¬ 
tual leg accommodates n copies of spin-l/2’s, and the 
symmetry class where xCq = XaCe = — 1 - While this 
is not the unique way to derive FQI, the imminent goal 
of this paper is to show how to produce an example of 
the featureless state for honeycomb spin-l/2’s, which is 
achievable with this particular choice of the symmetry 
class. One can choose the bond tensor to be the max¬ 
imally entangled state (g) icr2, where acts on 

the flavor space and 1(72 denotes the spin singlet formed 
by two virtual spin-l/ 2 ’s. For the site tensor, the most 
general form of a spin singlet (satisfying spin-rotation 
symmetry) and Kramers singlet (satisfying time-reversal 
symmetry) tensor is given by 

= E (^^7 (I + I ;;ta;/ 3 t 7 )) 

a,/3,7 

fo ^a/37 (I 1^5 ■ia-i/31^7) T | 1^a1^/3-i7)) 

fo ^a/37 (I 1^5 \ >l^at/3t7))^ • (4) 

Each element of C* is real to preserve the time reversal 
symmetry, and = 0 due to the SU(2) spin 

rotation symmetry. The first spin inside the ket sep¬ 


arated by the semicolon denotes the physical spin, the 
other three are virtual spins from each of the three legs 
for a given site, and a, / 3 , 7 label the flavor of virtual 
spins, with 1 < a, yd, 7 < n. For T) = 2 (a single virtual 
spin-1/2 per leg), there is no PEPS solution satisfying all 
of lattice symmetries, hence we turn to the simplest non¬ 
trivial case with D = A (two flavors of virtual spin-l/ 2 ’s, 
n = 2 ). Setting Wcr = os in Eq. ([ 3 ]), we find that, in 
order to meet the condition of invariance under the Cq 
and (7 symmetries, only the following two independent 
solutions for the site tensor are possible: 

id) = p(2| t;;2t42) - I t7it2;2) - I t742t2) 

+ 2| i;t24-it2) - 14-;tii2t2) - 14-;tit24-2)), 

id) =p(|t72ti;i)-|t7iti;2) 

+ U;t2;iti)-U;ti;it2)). (5) 

V stands for cyclic permutation of the virtual states. 
The general site tensor consistent with all symmetry 
requirements can be written as a linear combination 
+ C2A^‘^\ with arbitrary real coefficients 
Cl, C2. We claim that the topologically trivial symmet¬ 
ric PEPS state is obtained from contracting all virtual 
legs of the site tensors = ciA^^^ + 02A^^^ and bond 
tensors T^, for appropriate choices of (ci,C2). 

There are two special cases, ci = 0 or C2 = 0 , for which 
the PEPS wave functions have the emergent U(l) IGG. 
For T" = in Eq. each ket state has two out 
of the three virtual spins with flavor index equal to 2, 
and one virtual spin with the flavor index 1 . The state 
made from contracting A^^^ obviously preserves the flavor 
quantum number. A U(l) operation U{ 0 )^ defined by 
multiplying the (flavor index)=1 virtual spin by but 
not the (flavor index)=2 spin, changes the site tensor by 
the phase These gauge transformations form a U(l) 
group and result in low-energy fluctuations of U(l) gauge 
fields, which is known to be confining at long wavelengths 
in two dimenstions 0- The same argument shows that 
the state made out of A^‘^^ will likely describe the U(l) 
spin liquid. 

Based on extensive numerical analyses for general 
C2/C1 ratios, we conclude that a featureless state with¬ 
out topological order, spontaneous symmetry breaking, 
or emergent gauge symmetry has been found at around 
(ci,C2) = ( 0 . 9 , 0 . 1 ). First, topological entanglement en¬ 
tropy has been extracted by fitting the calculated entan¬ 
glement entropy for varying system sizes. As shown in 
Fig. [SJa) we find the extrapolated value - 0.05 consistent 
with the absence of topological order. To evaluate the en¬ 
tanglement entropy, we imposed the periodic boundary 
condition along the a2-direction of our ansatz wave func¬ 
tion and employed the boundary theory of PEPS [^.[^. 

Gorrelation functions were measured for spin, bond, 
vector spin chirality, and scalar spin chirality, in order to 
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(a) ribbon(Ly —>■ cxd) 


( 6 ) cylinder(La; —>■ oo) 



FIG. 2. (Color online) Schematic figure of the lattice geom¬ 
etry used to evaluate (a) correlators and (b) entanglement 
entropy. A green dot denotes one unit cell, formed by u (red) 
and V (blue) sites of the honeycomb lattice. 
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FIG. 3. (Color online) Numerical results on spin-1/2 PEPS 
for (a) Entanglement entropy as a 

function of Ly. Two linear fits are shown, based on numerical 
data at Ly G {3,4, 5,6} and Ly G {5,6}, respectively, (b)-(e) 
Plots of correlation functions {6 A = A — (A)) for (b) spin, (c) 
bond (Vi = Siu • Si^;), (d) vector chirality (xi = Siu x S^^;) and 
(e) scalar chirality (Ai = • Siu x S^^;) as a function of 

the distance Rij = |xi — Xj|/|ai |; x^ is the position vector of 
the z-th unit cell. System size is fixed at Lx = 37. 


determine the gapped nature of the state and to check 
for the absence of symmetry breaking. For convenience in 
numerical calculation, a unit cell was redefined so that u- 
and u-site tensors are directly connected within a unit cell 
such that Xi^u — ^i,v = (si2 — 2ai)/3, where is the 

position vector of zi(u)-tensor at the z-th unit cell. Expec¬ 
tation values of local operators {'ip\Oi\'ip) and correlators 
{ 2 p\OiOj\ 2 p) are obtained by employing the MPS-MPO 
compression method |24l426| to approximately contract a 
given tensor network. Two-site variational compression 
is adopted [iBl where the initial ansatz is given by the zip- 


up algorithm [27| . The reduced bond dimension d, which 
is equal to the number of singular values kept during 
the compression process of MPS, is such that the trun¬ 
cated singular value is less than 10“^. To minimize the 
numerical instability caused by lower and upper bound¬ 
aries of the tensor network wave function, compression is 
repeated until the convergence of MPS |4^n) — 0|^n-i) 
is achieved, where |4^n) is the nth stage MPS and O is 
the MPO. Number of unit cells in the ai direction is 37. 

As one can see in Fig.O all correlators decay expo¬ 
nentially in the fixed system size L^ = 37. We double- 
checked the exponential decay by varying L^ and fix¬ 
ing the position of operators Oi and Oj to be at Lx/3 
and 2Lx/3, respectively. Exponential decays are clearly 
observed for all correlators as a function of Lx in Fig. 
m Lattice and spin rotation symmetries are numerically 
confirmed as well. Finally, we extracted the correlation 
lengths as a function of reduced bond dimension d and 
found them to saturate for large enough d [1^ . Based on 
such overwhelming body of evidences, we conclude that 
our ansatz PEPS is a topologically trivial, fully symmet¬ 
ric and gapped quantum state or FQI. 
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EIG. 4. (Golor online) (a) Spin, (b) bond, (c) vector chirality, 
and (d) scalar chirality correlations as functions of at the 
separation distance equal to Lx/2. Exponential decays are 
observed in all correlators. 

We also measured the entanglement entropy for the 
U(l) states (ci = 0 or C2 = 0). In both cases, it depends 
on the number of flavor 1 states and flavor 2 states 
we put on the boundary virtual legs. Physically, we 
can interpret flavor states as U(l) gauge charges. For 
boundary conditions with different total flavor numbers, 
we end up with wave functions supporting different 
number of electric field lines along the length of cylinder. 
These states are orthogonal to each other, and in general 
give different entanglement entropies. 

Conclusion - We have identified an exemplary state of 
a featureless quantum insulator on the honeycomb lattice 
of spin-l/2’s, based on the methodical search scheme de¬ 
veloped in Ref. [H. Four distinct classes have been idem 
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tified as a result of our search. We propose a state whose 
physical properties are consistent with FQL Compared 
to previous works on FQI’s the present method 

offers a much more systematic way to classify tensor 
network states consistent with symmetry and topologi¬ 
cal constraints. The gapped liquid phase we constructed 
is intrinsically strongly interacting, as there is no way 
to adiabatically connect them to free electronic states. 
These results may thus be relevant for correlated elec¬ 
tronic materials on the honeycomb lattice. For instance, 
evidences for a putative spin liquid ground state have 
been reported for Ba3CuSb2 09 [28|,l29|, in which the spin- 
1 /2 Cu may form a honeycomb lattice . 

Acknowledgements - SHJ and YR are supported by the 
Alfred P. Sloan fellowship and National Science Founda¬ 
tion under Grant No. DMR-1151440. HYL is supported 
by the NRF grant (No.2015RlDlAlA01059296). CMJ is 
supported by the David and Lucile Packard foundation 
and National Science Foundation under Grant No. NSF 
PHYll- 25915. 


[1] X. Wen, Quantum Field Theory of Many-Body Sys¬ 
tems: From the Origin of Sound to an Ori¬ 

gin of Light and Eleetrons^ Oxford Graduate Texts 
(OUP Oxford, 2007), ISBN 9780199227259, URL 

https://books.google.com/books?id=lfxpPgAACAAJ, 

[2] E. Lieb, T. Schultz, and D. Mattis, in Condensed Matter 
Physics and Exactly Soluble Models (Springer, 2004), pp. 
543-601. 

[3] M. Oshikawa, Phys. Rev. Lett. 84, 1535 (2000), URL 
http://link.aps.org/doi/10.1103/PhysRevLett.84.1535 

[4] M. Hastings, EPL (Europhysics Letters) 70, 824 (2005). 

[5] C.-M. Jian and M. Zaletel, Phys. 

Rev. B 93, 035114 (2016), URL 

http://link.aps.org/doi/10.1103/PhysRevB.93.035114, 

[6] 1. Kimchi, S. Parameswaran, A. M. Turner, E. Wang, and 

A. Vishwanath, Proceedings of the National Academy of 
Sciences 110, 16378 (2013). 

[7] B. Ware, 1. Kimchi, S. A. Parameswaran, and 

B. Bauer, Phys. Rev. B 92, 195105 (2015), URL 

http://link.aps.org/doi/10.1103/PhysRevB.92.195105, 

[8] S. Jiang and Y. Ran, Phys. 

Rev. B 92, 104414 (2015), URL 

http://link.aps.org/doi/10.1103/PhysRevB.92.104414, 

[9] D. Perez-Garcia, M. Sanz, C. E. Gonzalez- 

Guillen, M. M. Wolf, and J. 1. Cirac, New 

Journal of Physics 12, 025010 (2010), URL 

http://stacks.iop.org/1367-2630/12/i=2/a=025010, 

[10] H. H. Zhao, Z. Y. Xie, Q. N. Chen, Z. C. Wei, J. W. 
Cai, and T. Xiang, Phys. Rev. B 81, 174411 (2010), URL 
http://link.aps.org/doi/10.1103/PhysRevB.81.174411, 

[11] S. Singh, R. N. C. Pfeifer, and G. Vi¬ 


dal, Phys. Rev. A 82, 050301 (2010), URL 

http://link.aps.org/doi/10.1103/PhysRevA.82.050301, 

[12] B. Bauer, P. Corboz, R. Oriis, and M. Troyer, 

Phys. Rev. B 83, 125106 (2011), URL 

http://link.aps.org/doi/10.1103/PhysRevB.83.125106, 

[13] S. Singh, R. N. C. Pfeifer, and G. Vi¬ 
dal, Phys. Rev. B 83, 115125 (2011), URL 

http://link.aps.org/doi/10.1103/PhysRevB.83.115125, 

[14] A. Weichselbaum, Annals of Physics 327, 

2972 (2012), ISSN 0003-4916, URL 

http://www.sciencedirect.com/science/article/pii/S000349161 

[15] S. Singh and G. Vidal, Phys. 

Rev. B 86, 195114 (2012), URL 

http://link.aps.org/doi/10.1103/PhysRevB.86.195114, 

[16] A. Kitaev and J. Preskill, Phys. 

Rev. Lett. 96, 110404 (2006), URL 

http://link.aps.org/doi/10.1103/PhysRevLett.96.110404, 

[17] M. Levin and X.-G. Wen, Phys. 

Rev. Lett. 96, 110405 (2006), URL 

http://link.aps.org/doi/10.1103/PhysRevLett.96.110405, 

[18] B. Swingle and X.-G. Wen, arXiv preprint 
arXiv:1001.4517 (2010). 

[19] N. Schuch, 1. Cirac, and D. Perez-Garcfa, Annals 

of Physics 325, 2153 (2010), ISSN 0003-4916, URL 

http://www.sciencedirect.com/science/article/pii/S000349161 

[20] P. Kim, H. Lee, S. Jiang, B. Ware, C.-M. Jian, M. Zaletel, 

J. Han, and Y. Ran, Supplementary Information (2016). 

[21] A. M. Polyakov, Gauge fields and strings^ vol. 140 (Har¬ 
wood academic publishers Chur, 1987). 

[22] J. 1. Cirac, D. Poilblanc, N. Schuch, and E. Ver- 

straete, Phys. Rev. B 83, 245134 (2011), URL 

http://link.aps.org/doi/10.1103/PhysRevB.83.245134, 

[23] T. B. Wahl, S. T. HaBler, H.-H. Tu, J. 1. Cirac, 
and N. Schuch, Phys. Rev. B 90, 115133 (2014), URL 
http://link.aps.org/doi/10.1103/PhysRevB.90.115133, 

[24] E. Verstraete and J. 1. Cirac, arXiv preprint cond- 
mat/0407066 (2004). 

[25] U. Schollwck, Annals of Physics 326, 96 (2011), 

ISSN 0003-4916, January 2011 Special Issue, URL 

http://www.sciencedirect.com/science/article/pii/S000349161 

[26] R. Oriis, Annals of Physics 349, 

117 (2014), ISSN 0003-4916, URL 

http://www.sciencedirect.coin/science/article/pii/S000349161 

[27] E. Stoudenmire and S. R. White, New Journal of Physics 
12, 055026 (2010). 

[28] K. V. Shanavas, Z. S. Popovic, and S. Sat- 

pathy, Phys. Rev. B 89, 085130 (2014), URL 

http://link.aps.org/doi/10.1103/PhysRevB.89.085130 

[29] A. Smerald and E. Mila, Phys. 

Rev. B 90, 094422 (2014), URL 

http://link.aps.org/doi/10.1103/PhysRevB.90.094422, 

[30] S. Nakatsuji, K. Kuga, K. Kimura, R. Satake, 

N. Katayama, E. Nishibori, H. Sawa, R. Ishii, M. Hagi- 
wara, E. Bridges, et al.. Science 336, 559 (2012), 
http://www.sciencemag.org/content /336/6081/559.full.pdf, 

URL http://www.sciencemag.org/content/336/6081/559. abstraci 



